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ON ; Abstract 
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We determine the restrictions imposed by primordial nucleosynthesis upon a heavy 

Oh! 

tau neutrino, in the presence of v T annihilations into Majorons, as expected in a wide 

■ 

class of particle physics models of neutrino mass. We determine the equivalent number 
of light neutrino species N eq as a function of m UT and the v T -v T -Majoron coupling 
5_i ■ g. We show that for theoretically plausible g values > 10~ 4 present nucleosynthesis 

observations can not rule out v T masses in the MeV range. Moreover, these models give 
N eq < 3 in the v T mass region 1-10 MeV, for very reasonable values of g > 3 x 10 . 
The evasion of the cosmological limits brings new interest to the improvement of the 
present laboratory limit on the v T mass which can be achieved at a tau-charm factory. 
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1 Introduction 



Despite great experimental efforts, the tau-neutrino still remains as the only one which can 
have mass in the MeV range. The present experimental limit on its mass is [[j]: 

m„ r <23MeV (1) 

Further progress will have to wait for the improvements expected at future tau-charm or B 
factories 0. On the other hand, many particle physics models of massive neutrinos lead 
to a tau neutrino with mass in the MeV range ||. Moreover such a neutrino may have 
interesting cosmological implications @]. It is therefore interesting to examine critically the 
cosmological constraints. 

The first comes from the critical density argument ||. However, as has been widely 
illustrated with many particle physics models where neutrinos acquire their mass by the 
spontaneous violation of a global lepton number symmetry |J, this limit can be avoided 
due to the existence of fast v T decays |7|, ||, § and/or annihilations JT0|, || into Majorons. 



Although the Majoron was first introduced in the context of the seesaw model |TIJ the 
spontaneous breaking of lepton number can be realized in many different models. There is 
only one important constraint on its properties following from the precision measurements 
of the invisible Z width at LEP, namely the Majoron must be mostly singlet under the 
SU(2) <S> U(l) symmetry. It has been noted that, in many models of this type the relic 
v T number density can be depleted well below the required value for all masses obeying eq. 

In order to demonstrate the cosmological viability of the MeV tau neutrino we must 



also consider the restrictions that follow from primordial nucleosynthesis considerations [|12| . 
In the standard model, these rule out v T masses in the range O, IT 



0.5 MeV < m VT < 35 MeV (2) 

This would imply that m VT < 0.5 MeV is the nucleosynthesis limit for the case of a Majorana 
tau neutrino. Here we will only assume that v T is a Majorana particle, which is the most likely 
possibility. This assumes for the maximum allowed effective number of extra neutrino species 
AN eq during nucleosynthesis either 0.4 or 0.6. Recent contradictory data on the primordial 
deuterium abundance [TJ, [UJ may cast some doubts on the validity of this assumption (for 
recent analysis see refs. [ 17 , |T8|]). In particular, if AN eq = 1 is allowed [|Tj|], there may be 



an open window for neutrino mass somewhere near 20 MeV. However it has been shown in 



ref. [19| that this window actually does not exist, when one carefully takes into account 
the influence of non-equilibrium electronic neutrinos on the neutron-to-proton ratio. These 
neutrinos would come from massive v T annihilations v T v T — > v e v e . 

However one knows that new interactions capable of depleting MeV v T density in the 
cosmic plasma are needed, at some level, in order to comply with the limit on the relic 
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Figure 1: Feynman diagrams for annihilations of tau neutrinos into Majorons. 



neutrino density. It is therefore reasonable to analyse their possible effect in relation with 
the primordial nucleosynthesis constraints |20| . 



In this paper we analyse the effect of neutrinos with large annihilation cross sections 
into Majorons. In order to compute the relevant annihilation rates we must parametrize 
the majoron interactions. These arise from the diagrams shown in Fig. (|l]). The t-channel 
diagram is present in all Majoron models, while the strength of the s-channel scalar exchange 
diagram is somewhat model-dependent. 

One way of writing the couplings of Majorons to neutrinos is using the fact that the 
Majorons are Nambu-Goldstone bosons and hence have derivative couplings. This is the so 
called polar coordinate method. The other method is to use a pseudoscalar interaction, some- 
times refered to as the cartesian method. The two methods are equivalent, even for second 
order processes as we are considering here, if we include all the Feynman diagrams contribut- 
ing at that order to the process of interest v T v T — *■ J J 0. In our calculations throughout this 
paper we will use the cartesian method of parametrizing the majoron interactions. Though 
we must in principle include also the s-channel diagram in Fig. (JJ), we will neglect this 
contribution. We explicitly show in the Appendix, that it is justified in our case to use only 
the t-channel contribution in order to derive a conservative limit on neutrino mass m UT and 
majoron coupling g. 

We have determined the restrictions imposed by primordial nucleosynthesis upon such 



§ Although equivalent, for models with a large number of scalars and where the Majoron is a linear 
combination of the imaginary parts of several fields, like the model of Ref. p3[ |, the cartesian method is 
more convenient. 



a heavy tau neutrino in the presence of v T annihilations into Majorons. We show that if the 
v T v T Majoron coupling constant exceeds g > or so, a large v T mass in the MeV range 
is allowed by the present upper bounds on the extra number of neutrino species. As a result 
one cannot rule out any values of the u T mass up the present laboratory limit of eq. ([!]). 0. 

We also show how such g values are theoretically plausible in the context of the most 
attractive elementary particle physics models where MeV tau neutrinos arise, and which are 
based upon the spontaneous violation of lepton number. 



2 Evolution of v T number density in the presence of v T 
annihilations 

Massive tau neutrinos certainly interact with leptons via the standard weak interactions, 
v T v T <-> v e ,^e^i e + e~, as assumed in refs. ||L3 , jl4| . Moreover, in many particle physics where 



neutrinos acquire mass from the spontaneous violation of a global lepton number symmetry 
H heavy neutrinos, such as the v T , annihilate to Majorons v T v T — > J J via the diagonal 
coupling 

C = i-gJiyJa 2 iy T + H.c. (3) 
where v T represents a two-component Majorana spinor, in the notation of ref. || |2T], 
This corresponds, in the usual four-component notation to 

£ = i-jgJv T ^v T (4) 

The corresponding elastic processes do not change particle densities, but as long as they are 
effective they maintain all species with the same temperature. 

We now comment on the cosmological bound provided by the critical density argument 
||. In order to be consistent with cosmological limits, the relic abundance of the heavy 
Majorana tau neutrinos must be suppressed over and above what is provided by the standard 
model charged and neutral current weak interactions, as well as those derived from Fig. ([]]). 
This happens automatically in many Majoron models, where neutrinos decay with lifetimes 
shorter than required by the critical density constraint || [7|, || [J. For example, in Majoron 
models of the seesaw-type a massive v T will typically decay with lifetimes shorter than 
the one required in order to obey the critical density bound, but longer than the relevant 
nucleosynthesis time, as illustrated in figure 18 of ref. ||. Another example is provided by 
the model of ref. j2^| . A v T lifetime estimate was given for this model in Fig. 1 of ref. . 



It is seen explicitly that a v T of mass in the MeV range of interest to us is expected to be 



^In fact, with a larger coupling constant g > 10 3 it may be possible for a stable MeV v T to obey the 
critical density limit, suggesting a possible role of v T as dark matter. 



stable on the nucleosynthesis time scale, but decays with lifetimes shorter than required by 
the critical density bound. This corresponds to a range of off-diagonal neutrino- majoron 
couplings 1CT 10 > g ff -diagonal > 1CT 13 , which naturally occurs in many models. 

For simplicity, we will assume from now in this paper that the massive zz/s decay with 
lifetimes shorter than required by the critical density bound, but are stable on the time 
scale relevant for nucleosynthesis considerations. The more general case where both decays 
and annihilations are simultaneously active on the nucleosynthesis time scale will be treated 
elsewhere P5] . 



2.1 Before Weak Decoupling 

Let us assume first that all species are interacting so that they have the same temperature. 
The evolution of the v T density can be found from the corresponding Boltzmann equation, 

h UT +3Hn UT = - £ (w) U - (nllfj^) (5) 

i=J,e,i> e , M V \ n i ) J 

In this expression (<Jiv) is the thermal average of the annihilation cross section times the 
v T relative velocity v. Using the convention for the momenta as in figure [l|, its value for the 
process v T v' T «-> x i x' i is f\ 

(<nv) = -4f^ / dU UT dU K dU Xt dU xl (2n) 4 5 4 (p + p' - k - k') 

x I M | 2 e - E ^ T e- E ^ T (6) 

Here we have assumed kinetic equilibrium amongst the different species, as well as Boltzmann 
statistics. By | M | 2 we denote the invariant amplitude obtained with the usual Feynman 
rules for Majorana neutrinos |9|, |21], p2[ , summed over all spins (and averaged over initial 
spins). Moreover we have set dU^ = d 3 pA/ (2ti) 3 2E Pa . 



Following reference |26| we express (<Jiv) as a single integral using the dimension-less 
variable x = m UT /T, 



M = 5 5 %72< S /°° 2 ds (s - AmUa^VSK,^) (7) 



where Ki(x) are the modified Bessel functions of order i (see for instance f27f ) and s = (p+p') 2 
is the invariant of the process v T v' T <-> Xix\. Using the new variable rj = 1 — 4m^/s instead 
of s, 

4x f 1 r] / 2x \ 

= Ki(xj Jo V J^W^ aM l( WT^ (8) 



HHere v = [{pp'f - mi T ) 1 / 2 / E p E p , . 



The cross-sections of the different annihilation processes are listed below. For annihilations 
to Majorons we have 

9 4 1 ~VU f 1 + Vv 



°M = — W i H—^ 1 ) - 2^1 • (9) 
JK " 128rrm 2 7? L VI - y /fjJ v 'J v ; 



128ir m 2 T r] <- v l — ■ Si /rj J 

where we have divided by 2! in order to account for identical Majorons in the final state 
and divided by 4 in order to account the v T spin factors. For the standard weak interaction- 
induced annihilations v T v T «-> fifi, in the limit of massless products we take 

**fo) = ^^(&L + &), (io) 

where b 2 L + b 2 R = 1/2 for i = u e<fl and b 2 L + b 2 R = 2((-l/2 + sin 2 6 W ) 2 + (sin 2 6 W ) 2 ) ~ 0.25 for 
z = e. 

One may write evolution equations analogous to eq. (§) for the other species present 
in the plasma, namely v e ^ and e^. However we assume that the weak and electromagnetic 
interactions are effective enough to keep u^'s and e's densities in their equilibrium values, 
rik = n e k q for k = u e;/x , e. Thus we are left with a system of just two coupled Boltzmann 
equations: 

< + 3Hn UT = - £ (a iV ) {u 2 Vt - «) 2 ) - (ajv) U - Kl) 2 -%) = S Ut (11) 

2 

hj + 3Hnj = (ajv) « - (n^) 2 ^) = Sj (12) 



7? 



Now let us briefly describe our calculations. First we normalized the number densities 
to the number density of a massless neutrino species, 77 — 0.181T 3 , introducing the quanti- 
ties r a = n a /n , where a = u T , J, and the corresponding equilibrium functions r e £. We then 
have for the time derivative of n a 

n a = f^- = S a - 3Hn a 

dn a dr a 3 
W = no df +ra T no 

or, equivalently, 

^=(--3ffr a )i-ir a (13) 
al v n ' T 1 

On the other hand, the time derivative of the temperature is obtained from covariant energy 
conservation law 

p = -3H(p + P) f = -3H{p + P) J ^ f (14) 



*The general formula is given in the Appendix, eq. 12 



where p is the total energy density and P is the pressure. Finally, as p = p(T,rj,r Ur ) we 
can rewrite 

dp dp dp drj dp dr Vr 
dT = dT + Wj'dT + dr~~dT ' 

and for the normalized particle densities one has 

dr Vr _ y (dp_ df^drj_ _df^dr^s _ 3_ 

dT ^ \dT + drj dT + dr VT dT > 1 ' 

drj ( dp dp drj dp dr Ur ^ 3 

W - ~^8T + dr-jlT + d^-&) ~ T rj (16) 
where, for a = u T , J, we have introduced 

Z*=^-;;(^-r a ) (17) 



p + P^3Hn 

The final Boltzmann system for the normalized particle densities is obtained from eq. (|15| ) 
and eq. (pl|) introducing the dimension-less variable x previously defined. Denoting r' = 
dr/dx, we have 

This system is valid as long as the tau neutrinos are coupled to the weak interactions. 
The following is the complete set of entries in equations (|I^) and ( 313 ) for the equilibrium 
quantities, total energy density and pressure, respectively: 

C = l -^x*h(x) , = - 

0.1817T 2 u ; ' J 3 

j_ j_ 37r %4 (+ , 1 . n nfi /2 ( x ) ^ 

P = Pu + Pe + P 7 + pj + pu T = —T 1 + -rj + 0.06a; 



10 V 6 ' V \ 

7T 2 4 / 1 I 3 (x) 



P = P„ + P e + P-, + P, + P„, = — T 4 ( 1 + -o + 0.062--^ r„ T ] . (20) 



In these expressions we have introduced the integral functions Ij, where j = 1,2,3, defined 
as 



h(x) = / du u 2 exp(— x\/ 1 + u 2 
Jo 

POO 

■^2 0*0 — / du u 2 Vl + v? exp(— xyl + v? 
Jo 



h(x)= du exp(-s Vl + u 2 ) (21) 

V 1 + u 2 



2.2 Past Weak Decoupling 



Once the u T 's decouple from the standard weak interactions, they remain in contact only 
with Majorons. Then one has two different plasmas, one formed by z/ r 's and J's and the 



other by the rest of particles, each one with its own temperature^ (denoted as T and T 7 ). 
Let us define now the variables 

_ «V _ "V- 

X ~ T ' V ~ T 7 

We assume that the photon temperature evolves in the usual way, y = Hy. The evolution 
equation of the v T and J number densities are now simplified versions of (|TTD and (|12l), 
because S Vt = —Sj, 

h VT + 3Hn Ur = -Sj 

hj + 3Hnj = Sj (22) 

or, in terms of r a 's, 

, Sj dy 

r = 

Ut n Hy dx 

r'j = ~r' UT (23) 

Due to the second equation, the Boltzmann system reduces to a single evolution equation say, 
for r Ur . However, one must still determine dy/dx which differs from unity because T 7^ T 7 . 
An equation relating y and x is obtained using the energy balance condition for the v T + J 
plasma. If p = p Ur + pj and P = P Vt + Pj, we can write 



p = -3H(p + P) , where H 



\ 3M ? 



(24) 



pi 



The expressions for r^ 7 , p and P given in equations eq. ( p0|) need to be modified in order 
to take into account the fact that there are two distinct temperatures T and T 7 . This leads 
to the following equation 



dy _ 


y 


7r 2 rj 
20 x 5 ' 
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fh(x) 

\h(x) 


(h(x)\\ ' 




dx 


3(0.06^Mrv T 

V h(x) 


+ 


n 2 rj_\ 
60 x J 


*(*-°-18£§) 



(25) 



Here we defined I±{x) = —dl2(x)/dx. 



In order to determine the final frozen density of v T which will be relevant during nucleo- 
synthesis we have to solve numerically the corresponding set of differential equations. Before 
weak decoupling these are ([18]) and (|T9| ) , while after decoupling one should combine eq. ( |23| ) 
and eq. (p5|) , with the initial conditions r a = r^, a = J,v T valid at high temperatures. 



In Fig. (fj) we show the results of our calculations of the asymptotic (frozen) values of 
r VT m UT as a function of m UT for the standard model (g = 0) and for the Majoron model with 
different g values. Note that in the standard g = case we agree with the previous results 
of ref. [14| but get somewhat larger values than those obtained in ref. [13|. We ascribe this 



small discrepancy to the use, in ref.|l3] of an approximate expression for the v T energies, 
rather than the exact ones. 



^Eventually the massless neutrinos will also decouple from the second plasma, while the e + e pairs will 
annihilate to photons, thus generating the well known T Va — T 7 difference. 




Figure 2: Frozen values of t Vt tti Vt as a function of m Vr for the standard model (g = 0) and 
for the Majoron model with different g values in units of 10~ 5 . 
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Figure 3: Effective number of massless neutrinos equivalent to the contribution of heavy 
v T 's with different values of g in units of 10~ 5 . For comparison, the dashed line corresponds 
to the standard model case when g = 0. 

3 Nucleosynthesis constraints on (rn UT ,g) 

In this section we use the results obtained for the v T number density in order to constrain 
its mass from nucleosynthesis arguments. The value of t Vt as a function of {m UT ,g) is used 
in order to estimate the variation of the total energy density p tot = Pr + p VT . In pR all 
relativistic species are taken into account, including Majorons and two massless neutrinos, 
whereas p Ur is the energy density of the massive z/ T 's. 

In order to compare with the standard model situation it is convenient for us to express 
the effect of the v T mass and that of the presence of the Majoron in terms of an effective 
number of massless neutrino species (N eq ) which we calculate for each frozen value of r T (m UT ). 
In reality, the true value of r T (m UT ) is always larger than its frozen value, and we have taken 
this into account in order to obtain reliable limits in the low v T mass region. 

In order to derive the nucleosynthesis limits, first we developed a simple code for the 
numerical calculation of the neutron fraction r n , as presented e.g. in ref. |[28|| , varying 
the value of N eq . Then we incorporated p to t to this numerical code and performed the 



integration of the neutron-proton kinetic equations for each pair of (m VT , g) values, where g 
is the coupling constant which determines the strength of the v T annihilation cross section. 
Comparing r n (m UT ,g) with r n (N eq ) at T 7 ~ 0.075 MeV (the moment when practically all 
neutrons are wound up in 4 He), we can relate (m UT ,g) to N eq . 



We repeated this calculation adapting Kawano's nucleosynthesis code [29 to the case 
of a massive tau neutrino, both in the standard model and the Majoron extension. We have 
found that both methods are in good agreement. The results for the numerical calculations of 
the equivalent number of massless neutrinos during nucleosynthesis with the use of Kawano's 
numerical code are shown in figure ^|. For comparison the case of g = is shown (dashed line). 
From figure ^ one can see that, in the asymptotic limit of very large m VT the annihilation 
into Majorons is very inefficient (see eq. (|9|)), so that the effective N eq value is larger than in 
the standard g = case precisely by a factor 4/7, which corresponds to the extra Majoron 
degree of freedom. Thus, if we take also g very large we get just N eq = 2 + 4/7 ~ 2.57. Of 
course this asymptotic limit is already experimentally ruled out by the Aleph v T mass limit 
IHII and thus is not displayed. For m Ur values in the range from 10 to 23 MeV or so, N eq can 
be made acceptable, provided g is raised sufficiently For the intermediate v T mass region, 
1-10 MeV, and g > 3 x 10~ 4 the model may even give N eq < 3, which is possibly supported 
by some of the observational data. 

Finally, in the small v T mass limit the energy density of v T is roughly the same as that 
of the massless v e or , so that all g values shown in the figure lead to the same asymptotic 
value N eq = 3 + 4/7 ~ 3.57, corresponding to the three massless neutrinos plus Majoron 
(instead of 2 + 4/7 for a very heavy v T ). However, it might be that observations eventually 
could lead to a tighter limit N™ ax < 3.57. In such event a simple way out is to have the 
Majoron out-of-equilibrium, which would require a very small g value, g < (2 — 3) x 10~ 5 , 
so that the production of Majorons through annihilations of v T 's would be negligible 
Should the observations eventually lead to an even tighter limit N™ ax < 3 the situation is 
qualitatively different, as it would raise a conflict with the standard model. A possible way 
to lower N eq below three provided by our model is to have a massive v T in the MeV range 
and with a relatively strong coupling with Majorons. Indeed, one can see from Figure |3] that, 
while it is not possible in the standard model to account for N™ ax < 3, it is quite natural 
in our model to have N™ ax < 3 for a wide range of intermediate tau neutrino masses and 
reasonable large values of the coupling constants g. 

In summary, one sees that all v T masses below 23 MeV are allowed by the nucleosyn- 
thesis condition N eq < N™ ax if N™ ax > 3.57, provided that the coupling between z/ r 's and 
J's exceeds a value of a few times 10~ 4 . This situation seems at the moment compatible 



with the experimental data, at least the He and Li determinations (18 



**Of course such m UT values are allowed by nucleosynthesis in the absence of v T annihilations. 




Figure 4: The values of g(m VT ) above each line would be allowed by nucleosynthesis if one 
adopts the N™ ax = 3, 3.4, 3.8, 4.2 (from top to bottom). 

It is instructive to express the above results in the m Vr — g plane, as shown in figure |J 
The region above each curve is allowed for the corresponding N™ ax . 



4 Significance of the Nucleosynthesis Limits 

There has been a variety of Majoron models proposed in the literature |J. They are attrac- 
tive extensions of the standard electroweak model where neutrinos acquire mass by virtue of 
the spontaneous violation of a global lepton number symmetry. Apart from their phenomeno- 
logical interest as extensions of the lepton and/or Higgs sectors of the standard model ||, 
Majoron models offer the possibility of loosening the cosmological limits on neutrino masses, 
either because neutrinos decay or because they annihilate to Majorons. The first and most 
obvious of these is the limit that follows from the cosmological density argument ]7|, Efl. As 
we saw in the previous section one can also place limits on a heavy tau neutrino with mass 
in the MeV range by using primordial element abundances. We have determined the restric- 
tions imposed by primordial nucleosynthesis upon a heavy tau neutrino, in the presence of 
v T v T annihilations into Majorons. Our results are completely general and may be compared 
to any bound characterized by an allowed value of N™ ax . Given any N™ ax value one can 
readily obtain the allowed regions of m UT and the Majoron coupling constant g as shown 
in Fig. (U). As an example, a recent model- independent likelihood analysis of big bang 
nucleosynthesis based on 4 He and 7 Li determinations has claimed an upper limit N eq < 4.0 
(at 95% C.L.) [0. From Fig. (f|) this would imply that all ra Vr masses are allowed, as 



long as g exceeds 10~ 4 or so. However we believe that, in the present state of affairs, one 
should probably not assign a statistical confidence to nucleosynthesis results, to the extent 
that these are still dominated by systematic, rather than statistical errors. Strictly speaking, 
what Fig. (ji]) really displays is the equivalent neutrino number N eq for various combinations 
of (m UT ,g) parameters that give the same helium abundance, rather than real limits. Of 
course, from these contours which contain the raw information an educated reader can judge 
which helium abundance should be considered plausible or not. 

We now illustrate in concrete models the fact that such values of the v T v T Majoron 
coupling > are theoretically plausible. Different models imply different expectations 
for the Majoron coupling constants g and the relation they bear with the v r mass m„ T . Our 



discussion so far is applicable to the simplest seesaw Majoron model of ref. |jn| . In this case 
one expects that M 



where mjj is a typical Dirac neutrino mass and Mr oc (a) is the Majorana mass of the 
right-handed SU(2) ® U(l) singlet neutrino. Clearly g values in the range required by 
nucleosynthesis are quite reasonable say, for mp ~ 1 — 100 GeV and Mr ~ 10 4 — 10 8 GeV. 
Moreover, it is a good approximation in this model to neglect the s-channel scalar exchange 
diagram of Fig. (|l|). 

There is a wide class of alternative Majoron models characterized by a low scale of 
lepton number violation P, |3(| ^TJ . These models are attractive because they lead to a wide 
variety of processes which may be experimentally accessible ||. In this case one expects 
a simple direct correlation between the mass of the neutrinos and the magnitude of the 
diagonal couplings of neutrinos to Majorons. The neutrino mass is simply the product 
of the Yukawa coupling g and the vacuum expectation value (a) which characterizes the 
spontaneous violation of the global lepton number symmetry ||, 

m = g(a) (27) 

From this it follows that for m Vr ~ 10 MeV and (a) ~ 100 GeV one obtains g ~ 10~ 4 . This 
situation is therefore characteristic of models where lepton number spontaneously breaks at 
the weak scale. 

There are more complicated models where the degree of correlation between the v T mass 
m Vr and the lepton number violation scale may be different and may involve more free pa- 
rameters. Just to give a concrete example of such models, let us consider the supersymmetric 
models with spontaneous violation of R parity [Q. These models lead to 

/ \2 

m ex -p— (28) 

where (a) is identified with the vacuum expectation value of the right-handed SU{2) ® U(l) 
singlet sneutrino and M SUSY denotes a typical neutralino mass. The expected values of 



(g,m VT ) are depicted in Fig. (|5D, obtained when one varies the other relevant free parameters 
over a theoretically reasonable range. 

For all models with low-scale lepton number violation we have shown, by doing the 
full calculation, that the overall annihilation cross section for v T v T annihilation into two 
Majorons can be enhanced by an order of magnitude with respect to our above simplified 
calculation which neglected the s-channel scalar exchange diagram in Fig. ([!]). Although 
this would allow us to weaken our limits, the effect on g would only be a factor 10 1//4 < 2, so 
that the limits derived in figure ^ could be relaxed by a factor < 2 in this class of models. 

As a last comment, we note that the limits obtained in our paper could also be tight- 
ened by including the influence of non-equilibrium electronic neutrinos (and anti-neutrinos) 
produced by v T v T annihilations on the neutron-to-proton ratio |19j but, again, the effect is 
quite small on the bounds derived on g. 

Last but not least, we must compare the limits obtained by primordial big bang nucleo- 
synthesis with those derived from astrophysics. A new light particle, like the Majoron, may 
have an important effect on stellar evolution and this allows one to place stringent limits 
on the strength of the interaction of such particles P2"| . In the case we consider here, the 
Majorons interact predominantly with a heavy v T (with the mass in MeV range), so its influ- 
ence may be noticeable in supernova explosions when the temperature reaches tens of MeV. 
The bounds on Majoron properties which can be deduced from supernova physics have been 
widely discussed |33| and recently analysed in ref. |32| (see also references therein). For 
example a Majoron with Yukawa coupling to electronic neutrinos in the range 10~ 6 — 10~ 3 
could be important for supernova physics. However in our model this coupling to v e is much 
smaller. A Majoron coupling constant to tau-neutrinos around 10 may be potentially in- 
teresting for supernova physics and will be discussed elsewhere. Here we only mention that g 



values larger than (a few)xlO _5 ym/MeV may be dangerous because the coupling is strong 
enough for abundant production of Majorons in high temperature regions in the supernova 
core and simultaneously small enough so that the mean free path of the produced Majorons 
is larger than the central stellar core. Still the coupling g > 10~ 4 seems to be allowed. 



5 Conclusions 

In this paper we have investigated the implications for primordial nucleosynthesis of a heavy 
tau neutrino in the MeV range, in the presence of sufficiently strong v T annihilations into 
Majorons. We have determined the effective neutrino number N eq , or equivalently the pri- 
mordial helium abundance, and studied the level of sensitivity that it exhibits when expressed 
in terms of the underlying v T mass m Vr and coupling parameter g, the relevant coupling con- 
stant determining the v T v T annihilation cross section. Given the fact that present nucleo- 
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Figure 5: Expected values of m Vr and g in model of ref. [23] 



synthesis discussions are still plagued by systematics, it is useful to interpret our results this 
way, rather than as an actual limit in the statistical sense. For each m Vr value, one can in 
principle identify the corresponding lower bounds on g for which the v T v T annihilations to 
Majorons are sufficiently efficient in order not to be in conflict with nucleosynthesis. More- 
over, in contrast to the standard model, these models can account for a value of N eq < 3 if 
the v T mass lies in the region 1-10 MeV, provided > 3 x 10 -4 . 

We have been conservative in determining the nucleosyhthesis limits, to the extent that 
we have neglected model-dependent contributions from s-channel Higgs boson exchange, 
given in Fig. ([!]). This seems reasonable from the point of view of the relevant particle 



physics models [11, pL 30, 31 



We have also concluded that, indeed, the required choice of parameters can be naturally 
realized in Majoron models both with weak and large-scale lepton number violation. As a 
result, for sufficiently large but plausible values of the u T v T Majoron coupling > 10~ 4 one 
can not rule out any values of the v T mass up the present laboratory limit based on the 
cosmological argument. This highlights the importance of further experimental efforts in 
laboratory searches for the v T mass. Improvements expected at a tau-charm factory are 
indeed necessary, since the primordial nucleosynthesis constraints on the v T mass can be 
easily relaxed in a large class of extensions of the standard electroweak model. 



Appendix 



Here we show why one can neglect the s-channel diagram of Fig. ([I]) in the determination 
of the nucleosynthesis bound on m VT and majoron coupling g. 

The total cross-section for the annihilation to Majorons that corresponds to s-channel 
and t-channel diagrams of Fig. ([!]) is given by 



64tt m 2 . 



e 2 v ^+(l-2e) 



1 — 7] 

2r] 



where the parameter e is defined by 
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Hi 



and g UU Hi, Qujj are the couplings relevant for the s-channel diagram of Fig. (|l|). In Eq. 
(j30l) the sum is over all the CP-even scalars present in the model. From its definition, one 
can see that e is proportional to the couplings vvHi and HiJJ. When e — ^ the s-channel 
becomes zero. 

The value of e depends very much on the model. For the pure-singlet majoron models 



with low lepton number violation scale considered in ref. |30[ there is a strict correlation 
between the neutrino mass and the lepton number violation scale. In this case one has e = 1. 
For seesaw models, with lepton number violated at a large mass scale, one has e < 1. For 
the supersymmetric model with spontaneous breaking of R parity [23J at the weak scale one 
can show that e typically lies in a range around the value 1/2. In our analysis we wanted to 
stay as much model independent as possible. In order to have an idea of the dependence of 
our results on e we define 



F(x,e)= f 
Jo 



1 



V 



v y/2 



ri)K x { 



2x 



(31) 



which is just the integrand of eq. (g) in section |2.1| . In Fig. (^) we plot the function 
F(x, e) for e = 0, 1/2 and 1. We see that the value e = represents a lower bound on that 
integral. For most models we would get a higher value. If we notice that the cross section 
is proportional to g A , that diference in F(x, e) would translate into a smaller value needed 
for g to satisfy the nucleosynthesis bounds. Therefore, one can obtain a model-independent 
and conservative bound by taking the worst possible case, which corresponds to e = 0. Due 
to the dependence of F on g A the bounds on g would not be too sensitive to the value of e 
in the range of interest. This justifies our simplified expression for aj used in eq. @. 
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Figure 6: The function F(x, e) for various values of e 
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